Experimental Realization of Nonadiabatic Holonomic Quantum Computation 
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Due to the geometric nature, holonomic quantum computation is fault-tolerant against certain 
types of control errors. Although proposed more than a decade ago, the experimental realization 
of holonomic quantum computation is still an open challenge. In this Letter, we report the first 
experimental demonstration of nonadiabatic holonomic quantum computation in a liquid NMR 
quantum information processor. Two non-commuting one-qubit holonomic gates, rotations about 
x- and z-axes, and the two-qubit holonomic CNOT gate are realized by evolving the work qubits 
and an ancillary qubit nonadiabatically. The successful realizations of these universal elementary 
gates in nonadiabatic holonomic quantum computation demonstrates the experimental feasibility of 
this quantum computing paradigm. 

PACS numbers: 03.67.Ac, 03.67.Lx, 03.65. Vf, 76.60.-k 



1. Introduction. Holonomic quantum computation 
(HQC) was first proposed by Zanardi and Rasetti [jj. 
In their original work, the twisting of eigenspaces of adi- 
abatically varying Hamiltonian was used to manipulate 
quantum states in a robust manner. Due to the geomet- 
ric nature, HQC is robust against certain types of control 
errors. Since control errors are one main obstacle to the 
realization of quantum computation, HQC has become 
one promising quantum computation paradigm and at- 
tracted increasing interests recently |2Hl8j. 

Early HQC is based on adiabatic evolution, in which 
states are encoded in degenerate eigenstates of a Hamilto- 
nian, and gates are accomplished by adiabatically varying 
the Hamiltonian along a loop in the parameter space. Be- 
cause of the adiabatic requirement, long run-time is natu- 
rally required in the parametric control in adiabatic HQC 
(AHQC). This not only limits the gate speed, but also ex- 
poses the system to environment for a long time, and con- 
sequently leads to decoherence and reduces the efficiency 
of AHQC. To overcome these drawbacks in AHQC, nona- 
diabatic HQC (NHQC) has been pursued, and several 
NHQC protocols have been proposed [l7l[i8|. In NHQC, 
the long run-time requirement is avoided, while still re- 
taining all the robust advantages, making NHQC a very 
appealing quantum computing paradigm. 

In this Letter, we report the first experimental realiza- 
tion of NHQC using a liquid NMR quantum information 
processor. The NHQC scheme we realize is based on a 
variant of the recently proposed NHQC scheme in Ref. 
18]. In our modified NHQC scheme, decoherence-free 
subspace is not used and nonadiabatic holonomic evolu- 
tion is achieved by nonadiabatically evolving the work 
qubits and an ancillary qubit. To experimentally realize 
universal quantum computation, nonadiabatic one-qubit 
holonomic rotation gates about x- and z-axes and the 
nonadiabatic two-qubit holonomic CNOT gate are suc- 
cessfully implemented using a three-qubit NMR quantum 
information processor. These results demonstrate the ex- 
perimental feasibility of NHQC. 

2. Theoretical protocol. We first briefly review the holo- 



nomic conditions. Consider an A^-dimensional quan- 
tum system with its Hamiltonian H(t). Assume the 
state of the system is initially in a Af-dimensional sub- 
space S(0) spanned by a set of orthonormal basis vectors 
{|<fe(0))}fcii- It has been proved that [U El the evo- 
lution operator is a holonomic matrix acting on 5(0) if 
\<f>k(t)) satisfy the following conditions, 
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= 0, k,l = l,...,M, 

period and \<f>k(t)) = 
Texp(— i J Q H(ti)dti)\4>k(0)) , T being time ordering. 

Now we construct the universal set of NHQC gates. 
For the nonadiabatic one-qubit holonomic rotation gates, 
a two-qubit system is used. We choose the logical qubit 
states as |0)l = |10), \1)l = |H}- By such a design, 
all the information of the logical qubit is encoded in the 
work qubit (the second qubit), and the first qubit acts 
as an ancillary qubit. We design two types of Hamil- 
tonians, H\(4ii) and H 2 ((j) 2 ), to respectively realize two 
non-commuting nonadiabatic one-qubit gates, 
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where a\ - 
J 2 sm(<£ 2 /2) 



Jicos(</>i/2), bi = Jisin(<fo/2), a 2 = 
b 2 = J 2 cos((j) 2 /2), I is one-qubit identity 
matrix, and X, Y, Z are Pauli matrices. In the ba- 
sis {| 00), 1 01), 1 10), 1 11)}, the evolution operators U^fa) 
and U 2 {t%) generated by H\{4>\) an d H 2 {<p 2 ) respec- 
tively read 
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where Ji7"i = and J2T2 = 7r. According to Eqs. 

©-(ini), it is readily to prove that both conditions (i) and 
(ii) are satisfied if the state of the two-qubit system is 
initially in the logical subspace 5f = {|0)l,|1)l}. So 
^? 1 ( t i) an( i U 2 2 (t 2 ) are holonomic matrices acting on 
5f. In the basis {|0) L ,|1) L }, f/f^n) and 0| 2 (r 2 ) are 
respectively equivalent to 



u xz {4>i) = 
u zx {4> 2 ) = 
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Then the one-qubit holonomic rotation gates about x- 
and z-axes acting on the space can be constructed by 
using U xz and U zx , 

R L Z {6) = U xz (0)U xz (~) -> C^(ri)^r f (n), (9) 

i?£(0) = E/«(0)tf, x (-^) -> C/ 2 (t 2 )Z7 2 -'(t 2 ). (10) 

From the above two gates, an arbitrary one-qubit NHQC 
operation can be built. 

The nontrivial two-qubit NHQC gate we realize is the 
nonadiabatic holonomic CNOT gate. A three-qubit sys- 
tem is used to implement this gate. 1 100) , 1 101), 1 1 10) 
and I 111) are encoded as |00) L , |01) L , |10) L and |11) L . 
We see that all the information of the logical two-qubit 
state is encoded in the two work qubits (the second qubit 
and the third qubit), and the first qubit acts as an ancil- 
lary qubit. The Hamiltonian H 3 for realizing the CNOT 
gate can be expressed as 



H3 =-^(Xi(7 2 — Z 2 )X 3 



X l {I 2 -Z 2 ){h 
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Letting the evolution time satisfy the condition 
J3T3 = tt/V2, the evolution operator in the basis 
{|000), |001), |010), |011) |100), 1 101) , |110), 1 111)1 reads 



0-3(73) =Diag[l, 1,1, -1,1,1, X] 



(12) 



According to Eqs. (fTTj) and (fT2|) . we can prove that 
both conditions (i) and (ii) are satisfied if the state 
of the three-qubit system is initially in the logical 
subspace 5 2 L = {|00) L , |01) L , |10) L , |11) L }. So U 3 (t 3 ) 
is a holonomic matrix acting on S 2 . In the basis 
{|00) L ,|01) L ,|10) L ,|11) L }, U 3 (t 3 ) is equivalent to the 
nonadiabatic holonomic CNOT gate. 

As the Hamiltonians H\(<f>{), H 2 (4> 2 ) and H 3 are time- 
independent, their holonomic evolution operators can be 



respectively written as 

Ut(r 1 )=U^ 1 ut (Ah), 
ut(r 2 )=Ul\ut(At 2 ), 
03(r 3 )=ng 1 3 (At 3 ), 
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where Ati (i £ {1,2,3}) is small time interval and its 
value is Ti/Ni, iVj being the number of the time steps 
of the holonomic evolution. By using a modification of 
the Trotter formula which is correct up to (At) 2 the 
short time evolutions respectively read 
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Here T-f^Ati), T$ 2 (At 2 ) and T 3 (At 3 ) can be realized by 
a combination of radio- frequency pulses and evolutions of 
the J-coupling constants between the neighboring qubits 
in NMR technique [20H26|. According to Eqs. (|T3> 
(TTBl . the nonadiabatic holonomic gates R^(9), R x (4>) and 
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Notably, in Eq. ([2"0|) . U 2 (t 2 ) can be implemented di- 
rectly with no approximations, thus it reads U 2 (t 2 ) = 

e i^L(X 1 I 2 -X 1 Z 2 ) _ 

3. Experimental procedures and results. Figures [T] (a) 
and (b) respectively illustrate the implementations of the 
nonadiabatic one-qubit and two-qubit holonomic gates. 
The diethyl-fluoromalonate dissolved in d6 acetone is 
used as the NMR quantum processor. 13 C, 19 F and 1 H 
nuclear spins respectively act as the ancillary qubit and 
the two work qubits. Both in the one-qubit and two-qubit 
cases, the ancillary qubit stays in state |1)(1| before and 
after the nonadiabatic holonomic evolutions. The input 
states (output states) of the work qubits are denoted as 
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FIG. 1: (color online) Circuits for the NHQC gates, (a) The 
one-qubit NHQC gates, (b) The two-qubit NHQC gate. Both 
in (a) and (b), 13 C acts as an ancillary qubit and stays in state 
|1)(1| before and after the nonadiabatic holonomic evolutions. 
19 F and 1 H nuclear spins are the two work qubits. 
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Input State 

FIG. 2: (color online) The experimental unattenuated output 
state fidelities for the NHQC gates. In (a), (b), (c) and (d) are 
the fidelities of p£ ut for R%(ir/2), -RJ(tt), -R£(tt/2) and R%(tt) 
respectively, applied to input states X, Y and Z. R^(n/2) 
and Rx(iv/2) are implemented on 19 F; R^(tv) and R^(tt) are 
implemented on X H. In (e) are the fidelities of p^ ut for U^ not , 
applied to 15 different input states listed in Eq. (|23[1 . The 
average fidelities (the red solid horizontal lines) are 97.6%, 
97.3%, 97.9%, 95.7% and 93.12% in (a)-(e), respectively. 
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3: (color online) Experimental spectra 


of 13 C. (a) 


and 



FIG 

(d) are respectively the spectra obtained by observing the 
states with X H and 19 F in pf n = X, with no holonomic oper- 
ations, (b), (c), (e) and (f) are the spectra of p^ ut , starting 
with the initial states pf n =X, and applying the holonomic 
operations R^(tt), R^(-k), R^(tv/2) and R^(n/2), respec- 
tively. Rx{k) and i?z(7r) are implemented on *H; R^(ty/2) 
and Rz(tt/2) are implemented on 19 F. (g) is the spectrum of 
the initial state pf n = IX. (h) is the spectrum of pf ut after 
the holonomic operation U^not with the initial state pf n = IX. 
(a), (d) and (g) are used as reference spectra, to which (b), 
(c), (e), (f) and (h) are normalized. All the observation is 
realized by transfering the states of the work qubits to 13 C 
and then observing 13 C. 



the ancillary qubit in Without loss of generality, 

the spectator work qubit is prepared in |0)(0|. Specifi- 
cally, the NMR processor is initialized in the pseudopure 
states |1)(1| ® pf n ® |0)(0| (for R%(tt/2) and R£(ir/2)) 
or |1)(1| ® |0)(0| ® pf n (for R%(tt) and R%(ir)). For the 
CNOT gate, the ancillary qubit is also prepared in |1)(1| 
and the whole state of the NMR processor is |1)(1| <g> pf n . 
In terms of the deviation matrices [30l | , the input states 
pf n and pf n are prepared in each of the following sets 

f&€{X,Y,Z}, (22) 



Pfn (Pout) and pf n (pf ut ) for the one-qubit and two-qubit 
cases respectively. 

We here realize the following four one-qubit NHQC 
gates, R%(ir/2), R%(tt), R%(n/2), R%(ir), and the NHQC 
CNOT gate U^ not . In order to demonstrate we can im- 
plement one-qubit NHQC gates on both 19 F and 1 H, 
Rz{tt/2) and R^(n/2) are implemented on 19 F and 
Rz(n) and iZ^W are implemented on 1 H. In our ex- 
periments, the number of iterations are chosen to be 
Ni = 3, N 2 = 2 and iV 3 = 2. We prepare the ini- 
tial states using the cat-state method [27H29j . For the 
one-qubit gates, we prepare the work qubit in pf n and 



pf n £{IX, IY, IZ, XI, XX, XY, XZ, YI, 

YX,YY,YZ,ZI,ZX,ZY,ZZ}. (23) 

The output states p^ ut and p^ ut are determined by 
quantum state tomography (QST) [3l|. To measure 
the sameness of the theoretical output state pth and 
the experimental output state Poyt, the attenuated and 
unattenuated state fidelities [32|, [33j] , which are respec- 
tively defined as Tr(p ou tPth)/ y / ^(p th p th )Ti(p in p m ) and 

Tr(PoutPth)/y/T^(PoutPout)T^(PthPth), are used. The at- 
tenuated fidelity takes into account the signal loss, while 



4 



the unattenuated fidelity ignores certain errors due to the 
signal loss and quantifies how similar in direction p out 
and pth are [H, |33| . The average experimental attenu- 
ated fidelities are 60.7%, 61.8%, 86.1%, 77.7%, 47.9% for 
the output states of R%(n/2), R%(n), R%(n/2), R%{n) 
and U^ not respectively, while their average experimental 
unattenuated fidelities are 97.6%, 97.3%, 97.9%, 95.7% 
and 93.12% respectively. These numbers are in-line 
with the results of other experiments done using diethyl- 



fluoromalonate 34]. The differences between the atten- 



uated and unattenuated fidelities are consistent with the 
signal loss rates measured in our experiments (see Sup- 
porting Material). Figure [2] shows the unattenuated out- 
put state fidelities in our expriments. Figure [3] shows 
example NMR experimental spectra. Figure [3] (a) (Fig. 
131(d)) is the 13 C spectrum of the input state pf n = X for 
*H ( 19 F), with 19 F ( X H) instate |0)(0|. Figures|3j(b), (c), 
(c), (f) show the 13 C spectra of the output states Pg Ut , af- 
ter implementing R^(n), R^(n), R£(n/2) and R^(tt/2) 
to the input states pf n = X, respectively. Figures [3] (g) , 
(h) show the spectra of pf n = IX and p^ ut after applying 
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FIG. 4: (color online) The QPT \ matrices of one-qubit holo- 
nomic gates ii«(ir), i?£(f ) and i?£(7r). The (a) and 

(c) columns are the real parts and imaginary parts of the the- 
oretical x matrices, respectively. The (b) and (d) columns 
are the real parts and imaginary parts of the experimental 
X matrices, respectively. The numbers in the x- and y- axes 
refer to the operators in the operator basis set {/, X, —iY, Z}. 



Quantum process tomography (QPT) [35J is used to 
quantitatively describe the implementation of the NHQC 
gates. According to QPT, each quantum process is char- 
acterized by a x matrix. For a given input state pi n , the 
output state is expressed as p out = 'Ek,iXki e kPin e l ■> where 
ef. belongs to an operation basis set. The elements of the 
operator basis set for the one-qubit and two-qubit cases 



can be respectively chosen as 

e k e{I,X,-iY,Z},k=l,...,4, 



(24) 



e k e {II, IX, -UY, IZ, XI, XX, -iXY, XZ, -iYI, 
- iYX, -YY, -iYZ, ZI, ZX, -iZY, ZZ}, 
fe=l,...,16. (25) 

The QPT x matrix is calculated using the output states 
via the technique described in Ref. |l9|. The experi- 
mental x f° r one-qubit and two-qubit gates are shown 
in Figs. 2] and [5] respectively. We use x fidelities to 
evaluate the performance of NHQC gates. The attenu- 
ated x fidelities |Tr(x e xpXtfc)l 03 > which take into ac- 
count the signal loss, are 70.5%, 71.3%, 89.5%, 83.3% 
and 51.2% for the R%(n/2), R^{n), R%(n/2), R^{n) 
and U^ not gates respectively. The deviations between 
Xth and Xexp are mainly caused by overall loss of sig- 
nal. To see the sameness of theoretical and experimental 
quantum processes when ignoring certain errors due to 
signal loss, we use the unattenuated x. fidelity defined 

as \Tr(xe Xp xl h )\/\/Tr(Xe XP xlp)T<Xthxl h ) 03 M 
The unattenuated experimental x fidelities of the gates 
^(tt/2), R%(tt), R%(n/2), RL(n) and U^ not are 95.9%, 
95.9%, 98.1%, 96.3% and 91.43%, respectively. It is in- 
teresting to note that the Trotter approximations in Eqs. 
([T^| , ([T7| and (fI5]l give very good approximations to the 
exact evolution and the theoretical x fidelities are 99.2%, 
98.6%, 99.2%, 97.4% and 98.7% for R^{n/2), R^(tt), 
R^(tt/2), RL(tt) and U^ not respectively. 




FIG. 5: (color online) The QPT x matrices of (7 c L „ ot . 
(a) and (c) are the real part and imaginary part of 
the theoretical x matrix, respectively. (b) and (d) are 
the real part and imaginary part of the experimental x 
matrix, respectively. The numbers 1 to 16 in the al- 
and y- axes refer to the operators in the operator ba- 
sis set {II, IX, -UY, IZ, XI, XX, -iXY, XZ, -iYI, -iYX, 
-YY, -iYZ, ZI, ZX, -iZY, ZZ}. 

4- Summary. As a proof of principle, we experimentally 
implemented NHQC via a NMR quantum information 
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processor using a variant version of the scheme proposed 
in Ref. (l8j . In our experiments, one-qubit nonadia- 
batic holonomic gates and two-qubit holonomic CNOT 
gate, which compose a universal set of NHQC gates, are 
implemented by using an ancillary qubit which provides 
the additional dimension needed in the holonomic evo- 
lution. This is the first experimental demonstration of 
NHQC, which is a step towards fault-tolerant quantum 
computing. The successful realizations of these universal 
elementary gates in NHQC demonstrate the feasibility 



of implementing NHQC using present experimental tech- 
niques. 
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